Dynamical mean-field approximation for pair contact process with a particle source 
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The one-dimensional pair contact process with a particle source is studied by using dynamical cluster 
mean-field approximations with sites up to n = 12. The results obtained for different levels of approximation 
become convergent especially for n > 6 and allow us to derive reliable extrapolations to the limit n ^ oo. 
At the zero source limit, the critical point exhibits a discontinuity whose magnitude vanishes with 1/n. The 
coherent anomaly analysis of data supports that the vanishing of order parameter and density of isolated 
particles has the same critical behavior. In contrast to an earlier prediction, the present approximation 
does not support the existence of critical behavior in the inactive phase where the frozen density of isolated 
particles depends on the initial state. 

PACS numbers: 02.70.Lq, 05.50.-|-q, 05.70.Ln, 64.60. Cn 
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The study of phase transitions from active fluctuating 
phases into absorbing states has attracted considerable 
interest in the last decade The static exponents of 

these phase transitions belong frequently to the directed 
percolation (DP) universality class, however, there are 
also examples for new universality classes Q . Open ques- 
tions are related to the necessary conditions that can de- 
stroy this robust universal behavior. One of the simplest 
model showing static DP behavior is the pair contact 
process (POP) which was proposed to realize a system 
involving infinitely many absorbing states [Q. Recently, 
Dickman et. al. introduced a modified PCP model to ex- 
plore the robustness of DP transition . In this model, 
each site of the one-dimensional lattice is either vacant or 
occupied by a single particle. A (randomly chosen) pair 
of nearest-neighbor particles is annihilated with a prob- 
ability p or an additional particle is created around the 
given pair with a probability 1 — p if it is not forbidden 
by double occupancy. In the extended model, an exter- 
nal particle source is introduced that attempts to insert 
isolated particles with a rate of h. This system exhibits 
an active phase when p is smaller than a critical value, 
Pc- For p > Pc the system evolves into a frozen (absorb- 
ing) state where the nearest- neighbor pairs are absent. 
Further details of the model can be found in Ref. [|. 
This extended model was studied by Monte Carlo (MC) 
simulations and dynamical cluster mean- field approxima- 
tion for quite large cluster sizes (ranged from n = 2 to 
6). Some disturbing behaviors, however, remained un- 
solved. For example, the analytical predictions tend not 
monotonously toward the MC results for h > when 
the cluster size is increased. Furthermore the analytical 
results shows a discontinuity in the variation of critical 
point Pc if /i — > 0. This observation is surprising because 
the present approximation has proved to be satisfactory 
in many cases for n < 6 |^ . It is expected that the further 
increase of cluster size will resolve these discrepancies. 

In this Brief Report, we discuss the results of dynam- 
ical mean-field approximations for cluster sizes as large 
as n = 12. The present approach allows us to give more 
accurate extrapolations to the limit n — > oo. Using this 



large number of data (from ri = 2 — 12) we can also im- 
prove the prediction of coherent anomaly method (CAM) 
[Q when considering the critical exponents around the 
transition point. Furthermore, we could study the frozen 
density of isolated particle concentration toward the sys- 
tem evolves in the inactive phase Q . 

To identify the transition point we have evaluated 
the p-dependence of the density of pairs of particles p2 
(henceforth considered as order parameter) as well as the 
density of isolated particles pi. The transition point is 
determined by the zero-point of p2 for h = and by the 
breaking point of the function pi{p) for h > 0. This latter 
criteria was particularly useful because of the extremely 
low value of p2 in a large interval of p. 

We shall not describe the details of the mean-field tech- 
nique because it has already been applied and demon- 
strated for the PCP model by several authors previously 
P,ppC|] . The derivation of the hierarchy of equations of 
motion for the configuration probabilities on the n-site 
clusters becomes complicated for large n. The technical 
difficulties increase drastically when enlarging the size 
of clusters. Traditionally, the set of equation of motion 
is solved numerically in the stationary states. Despite 
the slow convergency toward the stationary solution for 
large n, the numerical integration of the master equa- 
tions seems to be a more efficient method to find the 
stationary probability of configurations than the tradi- 
tional Newton-Raphson method. Using the numerical 
integration method we could determine all the configu- 
ration probabilities appearing on the 12-site cluster. (As 
an example, to get a data in last row of Table I requires 
4-week running on a personal computer.) It is worth 
mentioning that very recently this method is used suc- 
cessfuUy for the consideration of a stochastic sandpile 
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Our results are plotted in Figs. |l]- ||. Table I summa- 
rizes the predictions of all levels of approximations for pc 
and for the density of isolated particles at the criti- 
cal point at different values of h. Figure |l| demonstrates 
clearly that the predictions for pc become monotonously 
convergent for n > 6. The linear fit for n > 5 data 
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gives Pc = 0.075 in close agreement with the MC data 
iPc^'" = 0.077). The same good convergence can be ob- 
tained in the presence of external source. The compar- 
isons of linear fit with MC data are also listed in Table I 
(MC data are taken from Ref |]). 
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FIG. 1. Predictions for pc at different levels of approxima- 
tion (n = 2, . . . , 12) at /i = 0. Dashed line is the linear fit 
for n = 6, . . . , f 2 data and arrow shows the MC data. Insert: 
Discontinuity in pc at ft as a function of levels of approx- 
imation, n = ?>,... , f 2-point levels are plotted. Dashed line 
corresponds to a power law with an exponent oj « 1.666. 

Surprisingly, although the n-site approximations con- 
verge to MC data at arbitrary value of h, a jump can be 
observed in pc when considering the limit ft, — > for all 
71. The magnitude of this jump is defined as 



Ap" = limp'; (ft) 



p:{h = o) 



(1) 



The inset of Fig Q shows Ap" for different levels of the 
approximation. This log-log plot suggests that the jump 
decreases as a power low function, i.e., Ap" oc with 
an exponent lu « 1.666. The explanation of this value 
of uj remains to be clarified. This discontinuity of Ap" 
may also be observed for higher dimension versions of the 
model |l|]. 

A distinct improvement of convergency can also be ob- 
served for higher values of n if we consider the density of 
isolated particles (finat at Pc- The extrapolations, based 



on predictions of n > 5 data, are in excellent agreement 
with MC data. 

Such a large number of approximations makes possible 
to apply the CAM analysis introduced by Suzuki 0. In 
the vicinity of the critical point the order parameter (^2) 
and the density of isolated particles (pi) are estimated 
by p2 oc a^(pj? - p) and pi - 0„at oc a^p"^ - p), where 
pj? denotes the prediction for the critical point at n-point 
level. To estimate the /? critical exponents we have plot- 
ted the amplitudes of mean- field results (a" and Oj) as 
a function of 
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(Pc/P. 
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where Pc denotes the result of MC simulation. Figure || 
shows the mean-field amplitudes in the absence of exter- 
nal source (ft = 0) for different levels of approximation. 



10 



0.1 




0.1 



FIG. 2. CAM scaling of the critical mean-field coefficients 
for the order parameter (02, open squares) and the density 
of isolated particles (ai , crosses) . Dashed lines correspond to 
the P — 0.276 exponent. The values of a" are multiplied by 
1.1. 

Using different subsets of CAM-data the estimated ex- 
ponents are 0.221,0.246,0.276, and 0.312. The dashed 
lines correspond to the f3 — 0.276 DP exponent. Al- 
though the error bar is still large but the tendency to 
DP is straightforward. These results confirm the earlier 
MC observations ||]. Namely, the order parameter and 
the density of isolated particles can be described by the 
same power laws when the system approaches the critical 
point from the active phase. At the same time the CAM 
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analysis of mean-field data helps to understand why the 
exponents agree. Figure ^ demonstrates that these two 
quantities are proportional to each other in the vicinity 
of critical point. (In this plot 1.1a" is compared with a2). 

In the presence of source {h > 0) the earlier MC sim- 
ulation suggested a shghtly modified (3 — 0.287 expo- 
nent ||]. Our mean-field data become also convergent 
for n > 6 suggesting an exponent close to DP class 
(/3 = 0.274 ± 0.038 at ft = 10). However, the region in (5„ 
is so small that the error bar of the estimated exponent is 
four times larger than the predicted change of MC result. 
Therefore, this method is incapable to distinguish such 
close exponents. 
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FIG. 3. The stationary value of particle concentration in 
the inactive phase at 6-point level ai h — 0. The initial con- 
centrations for the different curves are indicated. 

As noted above, the large-cluster approximation makes 
possible to test the prediction of static critical behavior 
of particle concentration in the inactive phase . Such a 
critical behavior of a static quantity in the inactive phase 
has already been observed in a different model |jl^. In 
a previous MC study of PCP model it was suggested 
that the density of isolated particles follows a power law 
puat _ ^ _ p^^l^i for p > in the absence of ex- 
ternal source Following Marques et al. [||, we have 
calculated the absorbing state selected by the system's 
dynamics starting from the same homogeneous initial 
state. However, different initial particle concentrations 
were chosen to study the robustness of the final absorb- 
ing state in the inactive phase. The initial concentration 



is ranged from (/)„ to 0.99. The results of 6-point approx- 
imation are plotted in Fig. ^ for different values of initial 
concentration at /i = 0. It suggests that the absorbing 
state selected by the system's dynamics depends on the 
initial condition. Similar behavior may be observed for 
all levels of approximation. The observed deviation of 
curves increases further by increasing p and the level of 
approximation. Obviously, the stationary value of con- 
centration becomes independent of the initial condition 
in the active phase. 

We have tried to apply CAM analysis on the data ob- 
tained from the same initial concentration for different 
levels of approximation. Even in this case, the irregu- 
larity of mean-field coefficients at different levels of ap- 
proximation does not allow to extract critical exponent. 
Shortly, the present approximation does not support the 
existence of static critical behavior in the inactive phase. 
In the light of this prediction further intensive MC sim- 
ulations are suggested to clarify the possible existence of 
natural absorbing states. 

In summary, we have demonstrated that the present 
dynamical mean-field approximations (for large n) yield 
adequate extrapolation to the limit n ^ oo even in the 
presence of external particle source that causes nonana- 
lytical behavior in the limit /i — > 0. According to this 
approximation the discontinuity of pc displays power-law 
decay. The application of CAM analysis demonstrates 
that the densities of nearest-neighbor pairs and isolated 
particles are proportional to each other in the vicinity 
of the critical point and supports the same critical ex- 
ponents of these quantities. This behavior is certainly 
related to the fact that the extinct pairs can leave extra 
solitary particles behind. Although our analysis is proved 
to be a useful tool for investigating stationary states in 
the active phase it does not confirm the occurrence of 
an unambiguous critical behavior in the inactive phase 
where the composition of the frozen state depends on 
the initial state. We hope that the obscurity of station- 
ary inactive state will stimulate further MC simulations 
and theoretical investigations. 
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TABLE I. results of rt-site approximations 



n 


Pc{h = 0) 


Pc{h = 0.1) 


Pc(/l=10) 




(l){h^0.1) 


<j}{h = 10) 


2 


0.2 


0.66667 


0.6667 





0.5 


0.5 


3 


0.12774 


0.1820 


0.182 


0.23 


0.4615 


0.4615 


4 


0.11846 


0.1502 


0.1573 


0.21665 


0.5 


0.5 


5 


0.11757 


0.1440 


0.1579 


0.17989 


0.4580 


0.477 


6 


0.10272 


0.1166 


0.1353 


0.22393 


0.4383 


0.4482 


7 


0.10069 


0.1148 


0.1323 


0.21840 


0.4450 


0.4511 


8 


0.09692 


0.1094 


0.1265 


0.22121 


0.4390 


0.4476 


9 


0.09387 


0.1044 


0.1224 


0.22879 


0.4356 


0.4442 


10 


0.09218 


0.1027 


0.1198 


0.22838 


0.4351 


0.4438 


11 


0.09048 


0.1004 


0.1175 


0.23120 


0.4325 


0.4418 


12 


0.08925 


0.0989 


0.1156 


0.23269 


0.4316 


0.4416 


oo 


0.075 


0.080 


0.095 


0.242 


0.422 


0.433 


SIM 


0.077091(5) 


0.086272(15) 


0.09785 


0.241(1) 


0.421(1) 


0.433 
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